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O ^ Abstract 

f-n | In this article we study the phase transition phenomenon for the Ising 

model under the action of a non-uniform external magnetic field. We show 
that the Ising model on the hypercubic lattice with a summable magnetic 
field has a first-order phase transition and, for any positive (resp. negative) 
and bounded magnetic field, the model does not present the phase transition 
O^l ' phenomenon whenever liminf/ii > 0, where h = {hi) ie x d is the external 

k*" ' magnetic field. 

m 

^ ! 1 Introduction 

The Lee- Yang theorem |10] is one of the most revisited results in Statistical Me- 
chanics fT2| 113] . especially because of its application to the study of phase 
transition phenomena. One consequence of this theorem is that for any nonzero 
uniform magnetic field, i.e. h = (hi) ieZ d, hi = h € K\{0} for all i £ Z d and 
(3 = 1/kT, the ferromagnetic (J > 0) Ising model on Z rf has an unique Gibbs 
measure in the thermodynamic limit, independently of the boundary conditions. 

In this paper we consider more general models where the magnetic field h 
is not supposed to be uniform. For such models the Lee- Yang Theorem is still 
valid, and a natural question is to ask if for these models the Lee- Yang Theorem 
still implies the absence of phase transition. The question of the uniqueness of 
the Gibbs measure in a non-uniform positive magnetic field it was considered by 
Georgii [6] and Fontes and Neves [3]. They considered the model with a non- 
negative random field of positive mean, and proved that for all temperatures that 
there exists a unique Gibbs state. Here, in the Theorem 4 we prove uniqueness 
of the Gibbs states at all temperatures for all models for which lim inf hi > ( 
no average assumption is required for the fields hi), showing that the Lee- Yang 
Theorem, in this case, still implies the absence of phase transition. Although we 
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present a proof for Theorem 4 in the hypercubic lattice, the same argument works 
for any connected amenable quasi-transitive graph with bounded degree. 

It is well known that the Ising model with postive magnetic field can present 
phase transition depending on the graph structure of the model. Jonasson and 
Steif in |7] showed that the Lee- Yang Theorem does not imply the absence of a 
first-order phase transition for this model in any nonamenable graph. Basuev in pQ 
obtained the same result for a class of amenable but not quasi-transitive graphs. 
In this paper, assuming that the magnetic field decays to zero, we obtain the 
same result for hypercubic lattices, which are examples of amenable and quasi- 
transitive graphs. In other words, even with the magnetic field taking positive 
values at all sites of the hypercubic lattice, we prove that the Ising model can 
present a first-order phase transition. 



2 Preliminaries and Main Results 

Consider the distance between x and y on Z rf given by ||x — y\\ = Ylt=i \ x i ~ Vi\ 
and, for any finite A C denote by <9A the set of sites in Z d whose distance 
to A is equal to 1. The energy in A of each configuration a € 0, = {— 1,+1} Z 
satisfying the boundary condition uj € in A (uj b.c), that is, <7j = uj{ V i G Z d \A, 
is given by the Hamiltonian 

h a(<t) = ~\ ai(T i ~ hi(7i > w 

where denotes the set of ordered pairs in A U dA of nearest neighbors. For 
J > and h = (hi) i& %d the Hamiltonian defines a ferromagnetic Ising model with 
external field h. 

The Gibbs measure in A with uj b.c. is the probability measure on (fi, J 7 ), 
where T is the sigma-algebra generated by the cylinder sets given by 

/4 ,h '» = (2) 

if a satisfies the u b.c. and zero otherwise. The normalization factor is the 
standard partition function 

^=^e-WW (3) 



with the sum over all configurations a satisfying the uj b.c. 

We denote by Gp,h the set of Gibbs measures given by the closed convex hull 
of the set of weak limits: 

P,h,u = w _ lim Ijfw (4) 

where A n C A n+ i and uj runs over all boundary conditions. 

We say that there is a first-order phase transition when the set Gp,h contains 
more than one measure, see [2J [5]. In the case of a non-zero uniform field, the 
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Lee- Yang theorem can be used to prove that the analyticity of the pressure with 
respect to the parameter h is equivalent to saying that the system has no phase 
transition, see [TT] . 

We write h € ^ 1 (Z d ) when ||h||i = J2iez d 1^*1 < 00 an( ^ h e ^°°(Z d ) when h is 
bounded, that is, sup fcGZ d \hk\ < 00. 

The paper is organized as follows: Our main result is presented in Section 3, 
roughly speaking it states that, if the external field h is summable on Z d , that 
is, h G £ 1 (Wi d ), then the model has a first-order phase transition. This result is 
obtained by a Peierls-type argument, based on contours already used in \W\ \IE\ 
and a cancelation argument for clusters of contours with opposite signs. 

In Section 4, we adress the question of non-summable h. We show that, if we 
have lim inf hi > (similarly for lim sup hi < 0), then the system has no first-order 
phase transition. We remark that although this fact is pretty obvious, from the 
physical point of view, the arguments we found in the literature are based upon 
translation-invariance [9]. For a general £°°(Z d ) magnetic field it is not clear that 
we can use those techniques. By controlling the Radon-Nikodyn derivative we 
show how to deal with the non-invariant case. 

3 The expansion in contours and the phase transition 

The arguments in this section can be generalized for higher dimensions, we take 
dimension 2 for simplicity. We will describe roughly the contours in Z 2 , for de- 
tails the readers can see |16t I18j. This approach uses a bijection between finite 
collections of contours in Z 2 * = Z 2 + (1/2, 1/2) and configurations with some fixed 
boundary condition. Without loss of generality, we always suppose that the finite 
set A is a square. This fact helps us to assure that, if we have a finite number of 
contours in the dual set of A and already fixed a boundary condition, then there 
is a configuration associated to these contours, see lemma 2.1.2 [18] , 

Theorem 1. If the magnetic field h in the hamiltonian in (0) belongs to £i(Z 2 ), 
then the model presents a phase transition when J > 3||h||i. 

Proof: let A C Z 2 be a finite set and suppose initially that hi > V i £ Z 2 . It 
follows from the second Griffiths inequality that 

/ \P,h,+ ^ / \/3,0,+ , 0,0,+ / f -.in 
to) A > to/A =1 ~ 2 ^A ({ (J i = - 1 }), 

where is the null magnetic field. Using this lower bound and the following 
standard inequality (see p"> page 170 for details) 

4'°' + (to = -1}) < </3) := J>n + 3)3"- 1 e- 2 ^™, 

n=4 

since clearly c(/3) — > 0, uniformly in A when /3 — > 00, we get that 

lim to)^ h <+ = 1. (5) 
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The next step is to show that linig_ s . 00 (<7j) /3 ' h ' = —1. For this, it will be 
convenient to consider the Hamiltonian of the Ising model with - b.c. given by 

= " E ^i°3 - !) - E ^fa + ( 6 ) 

Note that this normalization does not change the measure //^' h '~. 

We will identify Z 2 with the subset of ]R 2 of integer coordinates. Fix a finite 
set A C Z 2 , to each site i of A we associate the dual plaquette p* (i) having i at its 
center. We call plaquettes the unit squares in M. 2 whose corners are in Z 2 *. The 
dual set A* of A is the subset of Z 2 * of the corners of p*(i), where i is some site of 
A. 

For each configuration a\ satisfying the - b.c. we will associate a family of 
contours as follows: for each pair of nearest neighbors sites in <ta, where we have 
opposite signs, we consider the unit segment e joining the two sites. The dual 
segment e* will be the unit segment orthogonal to e passing through the middle 
point of e, and joining the two sites of A* that are closest to that middle point. 

The union of these dual unit segments, called edges, will form closed curves in 
M? such that the sites in A* will have degree 0, 2 or 4 and, using some rule we can 
"cut" the corner when the degree is 4, see |18| for details. 

This process give us a bijection between the configurations satisfying - b.c. 
and the finite sets of compatible contours in A*, that is, closed self-avoiding and 
pairwise mutually avoiding contours in A*. 

Fixed the - b.c, for each contour 7 there is an unique configuration u 7 which 
has 7 as unique contour. We define the interior 0/7, int7, as the set of all i € Z 2 
such that, for <r 7 , we have a% = +1 and d(i,j) > 1, where d denote the Euclidean 
distance in K 2 . We will use the notation int7 for the set of all i£Z 2 for which 
o~i = +1 in <7 7 . The volume of 7 is the cardinality of int7, V0I7 = | int7|. 

Let a a be a configuration satisfying the - b.c. We have already seen that 
there is a finite set of contours associated to a a- For each contour 7 all spins at 
int7 \ int7 have the same value. We say that 7 is of type + (resp. type -), if the 
value of these spins are +1 ( resp. —1). 

We define over the set of signed contours a function £ given by 



£(7) 



exp -2/3J| 7 | - 2/3 £ hi , 7 of type 



i£ int7 



(7) 



exp -2/3J| 7 | + 2/3 £ K , 7 of type +, 



\ i€ int7 / 

where I7I denote the number of unit segments that compose the contour. 

Denoting by the partition function corresponding to the - b.c, it follows 
from a straightfoward computation that: 

^-=i+e^ e rW (8) 

n>l (7i,...,7„) k=l 

A* —compatibles 



Phase Transition in Ising Model with Non-Uniform Magnetic Fields 



5 



We remark that the value £(7) depends explicitly on the type of the contour, 
which is different from the usual cluster expansion. Since 

(#- = 2/,f-({ ff4 = +l})-l (9) 

if we prove that /%' h ' ({oi = +1}) — > 0, uniformly in A when (3 — > 00, we are 
done. 

Whenever <jj = +1 there exists a contour type +, denoted by 7, that involves 

the site i. We use the notation 7 i to indicate that 7 involves the site i. We say 
that 7 involves 7', and we write 707', if any site i £ Z rf involved by 7', it is also 
involved by 7. 

Before continuing the theorem's proof we need the following lemma: 

Lemma 2. Let {71,..., 7„} be a collection of A* -compatible signed contours. If 
h € ^(Z d ), tfien 

(n \ n / n \ 

JJ e -2/3J| 7fc | < -Q e(7fc ) < g^Hhll! "Q e -2/3J| 7fc | _ 
fc=l / k=l \k=l / 

Proof: For a fixed collection of A*-compatible signed contours {71 ,j n } the rela- 
tion involving determines naturally a partial order in this set. Let {7 n , . . . ,7r fc } 
be the set of all maximal elements with respect to this partial order. So we have 
that 

n k ( \ 

ne(7 fc )=n^) n ^ ), m 

where the product over an empty set is equal to one. If {j : 7 n 7^} = for all 
1 < I < k, the upper and lower bounds claimed in the lemma are straightforward 
since int7 n n int7 rm = whenever 1 < m < I < k. On the other hand, for 
each 1 < I < k such that the set { j : 7 n jj } is not empty, consider the graph 
Gi = (Vi,Ei), where the vertex set V/ = {7 n } U {7^ : 7 r , 7j} and the edges set 
E[ is the collection of all unordered pairs {7^- , 7; } such that if j m jj or 7 m 7; , 
then 7 m jj and j m 7;, i.e., there are no intermediate contours between 7^ and 

From the definition of the contours and the edges set it is easy to see that Gi 
is a rooted tree where the root is the most exterior contour, i.e. this contour is 
not in the interior of any other contour. For 7j,7 s € Vi, we denote by d{^j^ s ) 
the usual graph distance in G[ . Let TVi be the set of vertices in V\ that are in the 
last generation of the tree, i.e., the vertices with degree one on G\ and that the 
distance to the root is maximal. Consider the following partition of 

t(I) 

TV t = [j TV t (j), 
j'=i 

where t(l) is the cardinality of the set {7 G Vi : d(j,TVi) = 1} and for each 
1 < j < t(l), the sets TVi(j) = {'Jji, ■ ■ ■ ,1j rj (i)} are the maximal subsets of TVi 
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possesing the same parent 7p ; (j)- By rearranging, the product (fT0|) can be written 

as 



n 

1=1 



hi) 

n i e(7p J( o) n ^ i n 

i =1 V TeTViO') / 7G G i 



(11) 



where is the tree subgraph of Gi induced by the vertices 



mu{ 7 p l(i)) ..., 7 p t(1)(i) }). 



Observe that the value of the function £ at the root 7 n appears in the product 
n 7 eG; wnen the graph Gi is not empty Putting 



Si(j) = mtjp. {l ) \ \J int7, 



it follows from the definition of £, independent of the sign of the contour 7p-(i), 
that the following bounds hold: 



n 

U7fe:7feeTV ! (i)U{7p. C0 }} 



-2/3J|7fcl 



and 



^(7P j( 0) II £(7)<e +2 ^ 6S ^ 



\ 



n 

V {7fc:7fceTV I y)U{7 Pj . (/) }} 



-2/JJ|7fel 



Proceeding as above, for each Gi, we get new 5/(j)'s that are disjoint for the 
ones previously defined. So using the reasoning iteratively we finish the proof. The 
proof ends when each of Gi is only the root 7 n or the root r ) n and the elements 
connected to it. 



□ 

To finish the proof of the theorem we use the lemma above and, by the contour 
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representation, we get the following upper bounds for /j,^ h ' ({<7j = +1}): 
/4' h '-({^ = +l})</4' h '~«3 7 0.}) 



< 



£+ . i + En>i E, ;nl + ; a nLi to) 

\ A* —compatibles 

! + En>lE {71,-.7n} EK=1 £(7fc) 
A* —compatibles 

(^" h| " + a*. 7 . )n{7V , ^ h| " (re- 

\ A* —compatibles 

^ l|h|11 + E„>i E {71 ,..., 7 „ } e-mn, ( m=1 FwSi) 

A* —compatibles 

< ^^e 4 ^" 1 <^e X p(-2/5J| 7 | + 2/3 + 40||h||i J • (12) 



< 



70 

70i 

Using the last inequality and | 7 | > 4, we obtain 

AiA h '"({^ = +!}) < ^exp(-2/3(Jn-3||h|| 1 ))n3 n . (13) 

n>4 

Taking J such that J > 3||h||i, we conclude that ^' h '~({crj = +1}) — > when 
(3 — > oo, uniformly in A. 

The case h < (/tj < 0, V i G Z 2 ) with h G ^i(Z 2 ) follows from the first case 
and the FKG inequality. For an arbitrary field h in G t\(Jl?) we use FKG and 
reduce the problem to one of the previous cases. □ 

Corollary. For all the magnetic field h G £i(7i d ), the Ising model presents a 
firts-order phase transition. 

Proof: First replace the fields hi by zero, at all sites of a finite large enough 
region T, in order to have J > 3||h||i, where h denotes the modified magnetic 
field. Thus there are two different measures /i_ and by Theorem 1. Now, take 
the local function A(o~) = Eier ^i a i- F° r any local function / : £2 — > K we have 
that 

Analogously for the + boundary condition. 

Suppose by absurd that the model with the magnetic field h does not present 
a first-order phase transition, thus for any /3 > fixed and for all local function / 
we have (/)^' h '~ = (f)^' h ' + - Applying this equality for f(a) = e~^ A ^ a \ we get 
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So it follows from the above equalities that fi-(fe) = fi + (fe ) for any (3 > 
and for any local function /, which is in contradiction with Theorem 1 by taking 
(3 > sufficiently large and f(a) = Oie~ A ^ . □ 

Remark: In the above argument we show how we can work with contours in a 
non-symmetric set up and also illustrate that a finite-energy, (quasi-)local change 
in a infinite system does not change global properties, such as (non) uniqueness of 
the Gibbs state. In fact, the summability guarantees that this proof works but, 
it is not a necessary condition for the phenomenon. For the Semi-Infinite Ising 
Model on Z x Z + Frohlich and Pfister [3] showed that there are two Gibbs states 
at sufficiently low temperature with a constant magnetic field h only over the sites 
in the boundary of the lattice and, as we already mentioned before, Basuev [1] 
obtained a phase transition in a more interesting case, where the magnetic field is 
constant at all sites of the lattice Z 2 x Z_|_. 



4 Absence of Phase Transition 

In this section we show the absence of phase transition in the non-summable case 
under condition that the magnetic field satisfies lim inf hi > for positive fields 
and lim sup hi < for negative ones. Although the proofs of the Lemma 3 and 
Theorem 4 below are presented for the hypercubic lattice, they are the same for 
any quasi-transitive connected amenable graph with uniformly bounded degree 
due the Theorem 5 of [7]. For the proof of the Theorem 4, in this section, we 
follow close [3]. 

Lemma 3. There is no phase transition in ferromagnetic Ising models with uni- 
form nonzero external magnetic field outside of a finite volume. 

Proof: Let h the magnetic field and Ao C Z d a finite set such that, hi = h for 
all i £ Ao and h £ M.\{0}. To show absence of phase transition for this model it is 
enough to show that (cr^' 11 ~ = (cr i )' 3 ' h > + for all /3 > and i G Z d . 

Consider Ao = {k} C Ac Z d , from now ( . ^ h,UJ and Z^' h,u denote respectively 
the expected value and the partition function with respect to the Gibbs measure 
defined by the Hamiltonian (pQ) with boundary condition ui and constant magnetic 
field h. It follows from the definition that, for all i € A 

fa>f h '" = (a, • e^~ h ^)^^. (14) 

We know that the expected value (o~i ■ e^^ hk ~ h ^ C7k )^ h,UJ is independent of the 
boundary conditions in the thermodynamical limit. So, we need to show that 
Z^ h,U) /Z^'* 1 ^ it is also independent of the boundary conditions when A /*• "L d . 

In order to evaluate the limit of the above ratio we will need to consider 
boundary conditions in A\{/c}. We define u\ as (oj\)i = uji for all i € dA with 
(uj\)k = +1 and, Ct>2 by (u)2)i = for all i £ dA and {u2)k = — 1 so 

Z A _ e z a\m + e z a\« / 15 j 



Z A eP Z A\M + 6 A\{fc} 
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To show that the above expression in the thermodynamic limit is independent 
of the boundary condition oj, it is enough to show that Z^l ^? /^a\ {k} does no * 
depend on wi and 002- To do this we take i = k in (|14p and h k +00, and we 
obtain 



„«."(«7 l -+n) lim I M*-K) A ^> + A \w 

= +1}) I l + e^*-^ 



A A\{k} 



for all finite set A containing k. Then, 



h f ^i^J (16) 

A\{fcj 

and the argument follows by induction in |Ao|. □ 

Theorem 4. J/h 6 £ 00 (Z d ) is such that liminf igZ d foj > 0, i/ien i/ie Isinj model 
with external field h has no phase transition. 

Proof: We omit the parameter /3 since the argument is valid for all /3 > 0. 
Let be /i = supj gZ d/ij, h = liminf igZ d hi and e a positive number such that 
< e < h/2. We denote by A C Z d a finite subset containing {i £ Z d ; < e}, 
h\ the restriction of the external field h to the volume A, h\ the constant magnetic 
field in the set A with the value h, analogously for h and h £ ^ the constant magnetic 
field taking values h — e. It follows from Fundamental Theorem of Calculus that 
for any finite T C 1> d with A C f, we have that the difference 



is equal to 



P {o-i;o-j) r dx + (ai) r - {<ji) T 

Jh-e .rrr A 



- ier\A 

By lemma [3] we know that there is no phase transition for the ferromagnetic 
models with uniform nonzero magnetic field outside of a finite volume. Then, by 
the FKG inequality taking the limit when T /Z d we have: 

< ((7i }^\A'+ - <^ a '^\a- = lim /if V (a t ;a 3 ) h r A ^ + dx 



Now, taking the limit when A Z d we obtain: 



£ jer\A 
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< (<7i) h ' + - (ai) h >~ < lim lim p [ V fa; <r^ A ' XnA ' + 
Since the truncated correlation functions are non-increasing in hk(k 6 Z d ), we 

get 

< fa) h '+ - fa) h '~ < 



< 



The last inequality comes again from FKG. By standard arguments [8], there 
is only one Gibbs state for the model. The analogous result holds when the mag- 
netic field is negative. | — | 
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lim lim /3 

A/'Z d r/ , Z d Jh-e 



\Gi](jj) r ax 



■ ier\A 

^A>r\A:+ 



lim lim (<Tj) r 
a />z d r /•i d 

lim (o-i)x A ' + - (o-j) 




h—e, 
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